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Numerical modeling of correlated single-electron tunneling 
in uniform 2D arrays of small conducting islands separated by 
tunnel junctions shows the possiblility of soliton-antisoliton 
avalanches. Though the time duration of any avalanche, and 
the total charge AQ = ne transferred across the array during 
the avalanche, are always finite, in arrays with length A'^ larger 
than certain critical value Nc and large width M S> f{N), the 
avalanche magnitude n may be exponentially large, resulting 
in particular in a giant increase of shot noise. Thermal fluc- 
tuations and disorder gradually suppress the avalanche effect. 



PACS numbers: 72.70-fni, 73.23.Hk 



During the past few years, there was much interest in 
correlated single-electron transport |l|,y] in large arrays 
of small conducting islands separated by tunnel barriers - 
see, e.g., Refs. P-|2^ and references therein. This atten- 
tion is due to several interesting properties of the arrays, 
the most basic of which is the existence of single-electron 
solitons ^,^. Such a soliton consists of an additional 
single-electron charge placed on an island of the array, 
surrounded by a group of totally neutral but strongly 
polarized islands which screen the electron field at large 
distances. Similarly, removal of an electron from a single 
island creates an "antisoliton" (i.e., a single- hole soliton) . 

The concept of single-electron solitons and antisolitons 
and their interaction allows a natural explanation of all 
the peculiarities of the array statistics and dynamics, in- 
cluding the Kosteritz-Thouless-like phase transition in 
the 2D case [p|J7[l^, single-electron oscillations of fre- 
quency / = I/e |i,|,13|. Coulomb drag HJl^, effects of 
disorder |12 l^,n8£0], and shot noise suppression pl|,p2[. 



The goal of this Letter is to report the prediction of a 
new phenomenon in 2D SET arrays: soliton/antisoliton 
avalanches which lead to shot noise enhancement rather 
than suppression. 

We have obtained theoretical evidence of this effect 
during Monte Carlo simulation of 2D array dynamics 
within the "orthodox" model of single-electron transport 
y . In this model, which is quantitatively valid when the 
tunnel conductances G between the array islands are suf- 
ficiently low (G ^ e^/h), single electron tunneling events 
are treated as incoherent transitions with rates 

r = ^ ^ (1) 

e2 1 - exp(-AM^/fcBr) ' ^ ' 

where AW is the drop of the electrostatic free energy 



caused by the particular transition. The energy W of 
a charge configuration was calculated within the usual 
approximation MM (strictly correct for an array suf- 
ficiently close to a conducting ground plane) in which 
the capacitance matrix includes only diagonal terms Cq 
and near-diagonal terms C which represent island self- 
capacitances and mutual capacitances between the neigh- 
boring islands, respectively. We studied not only the av- 
erage current / through the array, but also the spectral 
density Si{uj) of current fluctuations. 

The results show that at larger currents the Fano factor 
F = Si{0)/2el decreases to 1/A^ {N is the longtitudinal 
number of junctions), just as in ID arrays |21|, so that 
the low- frequency shot noise is indeed suppressed in com- 
parison with the Schottky value 2el. However, rather un- 
expectedly, we have found that in uniform arrays at low 
temperature T and low currents / (near the Coulomb 
blockade threshold) the Fano factor may be much larger 
than one - see, e.g., Fig. 1. 
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FIG. 1. Fano factor vs. dc current for several rectangular 
arrays of (A'^ — 1) x M similar islands, with the same length 
iV = 5 but various widths M, for Co = and T = 0. Point 
size corresponds to the Monte Carlo simulation accuracy; lines 
are only the guides for the eye. 

This fact invites the suggestion that electrons are 
transferred through the array in "bundles" with total 
charge AQ ^ e; however, since single-electron solitons 
repel P,H , the mechanism of the bundle formation was 
not immediately clear. The solution to this puzzle has 
turned out to be very simple. It is illustrated by Fig. 2 
which shows a few snapshots of a typical Monte-Carlo- 



simulated process at zero temperature, when the dc volt- 
age V applied to the array is just slightly (in this par- 
ticular case, by lO^^e/C ) above the Coulomb blockade 
threshold Vt- AtVfvVt, the soliton entrance into the ini- 
tially empty array is the main transport sequence bottle- 
neck, and takes place only after a considerable time pause 
when the array is free of any extra charges. After the en- 
trance from one electrode, the soliton starts to drift, hop- 
ping along the electric field applied to the array (Fig. 2a). 
When the soliton approaches the opposite electrode, its 
field may induce the entrance of one or more antisolitons, 
i.e. single-hole solitons (Fig. 2b). After these antisolitons 
have passed their transport bottleneck at the entrance, 
they drift mostly along the applied electric field (in the 
direction opposite to that of the soliton), the attracting 
field of the initial soliton having little effect on the drift. 
Because of this, soliton-antisoliton pairs frequently miss 
their chance to recombine, though lateral tunnel junc- 
tions between island rows allow such recombination and 
we do observe such events in our simulation. In turn, each 
antisoliton approaching the array electrode may trigger 
the entrance of one or more solitons (Fig. 2c), etc. This 
chain reaction results in a soliton-antisoliton avalanche, 
very much similar to an electric discharge in a gas due 
to surface impact ionization |23|. Notice that in ID ar- 
rays the soliton and antisoliton cannot pass each other, 
so that the recombination always happens; as a result 
avalanches are absent and F < 1 [pl[. 
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FIG. 2. A typical start of an avalanche in a 15 x 20-island 
array biased by dc voltage V ~ Vt — 2.35e/C (The voltage is 
applied between the array electrodes shown as hatched rect- 
angles.) Black squares mark the islands with an extra electron 
(the single-electron soliton centers), and white squares those 
with an extra hole (the antisoliton centers). Arrows show the 
directions of soliton and antisoliton motion. Co/C = 0.01, 
T = 0. (a) t = 8.48, (b) t = 30.6, (c) t = 119; here t is time 
measured (in the units of C/G) from the (random) moment 
of the first soliton entrance. 

Due to the stochastic character of transport within 
the framework of the orthodox theory, there is always 
a chance that every soliton and antisoliton leaves the ar- 
ray without triggering the entrance of any solitons of the 
opposite charge. Because of this, the avalanches always 
have a finite duration in time and finite "magnitude" n 
(the latter may be defined by the equation AQ — ne for 
the total charge transferred through the array during one 



avalanche). Nevertheless, the avalanches may be rather 
large: in wide arrays we have observed n's up to 4 x 10®, 
limited only by the computer simulation time. 

In the range where avalanches are distinct, computa- 
tion of F may be sped up considerably by using formulas 
for "charge blocks" , which were derived in Ref . p4] (in a 
different context): 

Sm = A (e2(„2) + j2^^2^ _ 2,j(„^)) ^ (2) 

I^e{n}/{r). (3) 

If avalanches do not overlap, n is just the magnitude 
of a single avalanche, while t is the interval between the 
beginnings of the adjacent avalanches. In the particular 
case / ^ 0, Eqs. @ and (|) yield 



F = 



in) 



(4) 



Figures 1 and 3 show the dependence of F on the main 
parameters of the system: array length N, width Af , is- 
land capacitance ratio Cq/C, and applied voltage V [the 
last dependence is presented parametrically via the in- 
duced dc current /(V^)]. Most of these dependences may 
be readily understood, at least qualitatively, using the 
avalanche picture discussed above. For example, if the ra- 
tio Co/C is increased, the soliton radius decreases ||,g] , so 
that the soliton-antisoliton interaction is quenched, and 
F is decreased. Larger N gives the soliton more time to 
induce solitons of the opposite sign, so that F grows. 
Finally, in very narrow arrays (small M) the soliton- 
antisoliton recombination suppresses the avalanche mag- 
nitude, so that the Fano factor also decreases. 
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FIG. 3. Fano factor at / — > (V -^ Vt) as a function 
of the array width, for T = and several values of N and 
Co/C. Open points: Monte Carlo simulation of the initial 
problem. Closed triangles: Monte Carlo simulation of the 
"macromodel" . Inset: Parameter p of the "macromodel" as a 
function of A^ and Co/C. 



We have found that the avalanche statistics depen- 
dence on M in the opposite hmit of large M may be 
understood as follows. In wide arrays, the probability p 
that a soliton triggers the entrance of an antisoliton, and 
does not recombine with it, should not depend much on 
the array width M . The average number / of rows sep- 
arating the soliton and the neighboring antisoliton, may 
depend on TV and Cq/C, but should be also virtually 
independent of M. This is why we may introduce an ap- 
proximate "macromodel" which is schematically shown 
in the inset of Fig. 4: an array of width M is presented 
as a parallel connection of m = M / f channels of equal 
width / = /{N^Cq/C). We break the avalanche into 
time stages of equal duration and assume that a passage 
of a soliton in some channel at stage i triggers the anti- 
soliton passage in each of the two neighboring channels 
at stage {i + 1) with probability p = p{N,Co/C). The 
lateral (open) sides of the array are described by setting 
the corresponding probability to 0. Mathematically, our 
macromodel is exactly the problem of the directed bond 
percolation problem on a 2D square lattice |^ within a 
stripe of width m. 
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FIG. 4. Avalanche magnitude hystogram for several arrays 
of different widths M, all with TV = 5 and Co ~ 0. Points: 
Monte Carlo simulation of the initial problem. Solid lines: 
Monte Carlo simulation of the "macromodel" (shown in the 
inset) for p = 0.512 and / — 3. 

Despite the approximate character of the macromodel, 
it allows the results of the initial problem to be re- 
produced remarkably well. For example. Fig. 4 shows 
avalanche amplitude histograms (probability to find an 
avalanche with particular magnitude n) for 5 arrays of 
different widths. The Monte Carlo simulation results 
virtually coincide with the macromodel results, despite 
only two parameters {p and /) were available for fitting 
all 5 curves. (An equally good fit was obtained for all 
other values of N and Cq/C we have studied if M was 
large enough, m ^ 1.) At the same time, the macro- 



model is much faster for simulation, so that with the 
same computer resources, results may be obtained in a 
much broader range of array parameters. For example, 
open points in Fig. 3 show Monte Carlo results for the 
initial model; for N ~ 10, we could hardly get accept- 
able accuracy at M > 40. For the macromodel described 
above, calculation of i^ in that point (with the same error 
bars) took 35 times less CPU time, and we could continue 
calculations all the way up to Af — 55 (closed triangles). 

The macromodel also gives a clear explanation why the 
growth of the avalanche magnitude with array width M 
saturates for shorter arrays, but is unbound for longer 
arrays (Fig. 3). The directed bond percolation problem 
on square lattice has a percolation threshold Pc = 0.6447, 
beyond which there is a finite probability of having an in- 
finite percolation cluster on an infinite lattice. We have 
found that for our problem, p becomes larger than pc if 
the array length is above a certain value Nc- {Nc « 8 
for Co = 0, and grows with the increase of the Cq/C 
ratio; for large enough Cq/C, p never reaches pc at all 
- see inset in Fig. 3.) This means that at TV > Ac and 
unlimited m there is a finite probability of having an infi- 
nite avalanche which, once started, would never end. The 
large, but finite, array width (and hence a large but finite 
m ^ 1) stops the avalanche growth in time, and limits 
its magnitude at a finite (though exponentially large in 
m) level. 

The avalanche effect is most strongly expressed in uni- 
form arrays at zero temperature, while both thermal fluc- 
tuations and array disorder gradually suppress it. Fig- 
ure 5 shows a typical dependence of the Fano factor on 
the applied voltage for several values of temperature. For 
T = 0, transport is possible only above the Coulomb 
blockade threshold V > Vt, so that below this point F 
is undetermined. However, even very small temperature 
fluctuations reveal the second branch of this dependence, 
since the initial energy barrier for soliton entrance to the 
array may be overcome by thermal activation, so that dc 
current becomes finite (though may be very small) . Fig- 
ure 5 shows that for V < Vt, the Fano factor decreases 
as applied voltage is decreased, approaching 1 for smaller 
V. The reason for this decrease is that as F ^ 0, it is 
harder and harder for the initial soliton to trigger the 
antisoliton entrance, so that single-soliton passages be- 
come the dominant component of transport. (At very 
small voltages, eV/2 < ksT, the Fano factor starts to 
grow again as 2kBT/eV due to quasi-equilibrium thermal 
fiuctuations. However, in our case of low temperatures, 
ksT <^ eVt, this growth corresponds to exponentially 
low dc current.) Larger temperature leads to gradual de- 
correlation of the moments when solitons and antisolitons 
enter the array, and hence to a gradual suppression of the 
avalanches which depend on this correlation. 
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FIG. 5. Fano factor as a function of applied voltage V for 
an array with A*' = 5, M = 10 and Co = 0, for several temper- 
atures (measured in e /ksC). Vt = 0.883e/C is the Coulomb 
blockade threshold. 

The effect of disorder is qualitatively similar (in qual- 
itative agreement vi^ith the results for the directed bond 
percolation p9]). For example, in a 5 x 20 junction ar- 
ray with Cq — 0, the introduction of background charges 
Qo — 0.5e to just 6 percent of the islands at random 
usually suppresses avalanches completely, reducing F to 
1. 

To summarize, we have found strong theoretical evi- 
dence for (and a simple explanation of) a new effect of 
single-electron-soliton and antisoliton avalanches in 2D 
arrays of small conducting islands. This effect arises due 
to soliton-assisted entrance of antisolitons into neighbor- 
ing, but different rows of the array, thus avoiding re- 
combination. The basic physics of the effect is well cap- 
tured by a simple "macromodel" , equivalent to the di- 
rected percolation model on a square-lattice strip of fi- 
nite width. Thermal fluctuations and array disorder sup- 
press the avalanches. Presently, we are working toward 
a quantitative characterization of the latter suppression, 
in particular with the goal to evaluate feasibility of ex- 
perimental observation of the soliton avalanches using 
various techniques of array fabrication [|5|,[7| p^|Jl5|Jl^ , |2Q | . 
Apparently, this may be easier to accomplish in an array 
with tunnel- uncoupled rows (essentially, a generalization 
of the 2-row system studied earlier Pjlq] ) , where tuning 
of dc voltages applied to each row may be sufficient to 
compensate for a moderate disorder. 
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